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Abstract 

The response functions for small spatial perturbations of a homogeneous granular fluid have been 
described recently. In appropriate dimensionless variables, they have the form of stationary state 
time correlation functions. Here, these functions are expressed in terms of reduced single particle 
functions that are expected to obey a linear kinetic equation. The functional assumption required 
for such a kinetic equation, and a Markov approximation for its implementation are discussed. 
If, in addition, static velocity correlations are neglected, a granular fluid version of the linearized 
Enskog kinetic theory is obtained. The derivation makes no a priori limitation on the density, space 
and time scale, nor degree of inelasticity. As an illustration, recently derived Helfand and Green- 
Kubo expressions for the Navier-Stokes order transport coefficients are evaluated with this kinetic 
theory. The results are in agreement with those obtained from the Chapman-Enskog solution to 
the nonlinear Enskog kinetic equation. 

PACS numbers: 05.20.Dd,45.70.-n,05.60.-k,47.10.ab 
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I. INTRODUCTION 



One of the most productive methods to study transport in normal fluids is through the 
measurement, simulation, and theory of linear response functions [1-4]. Of particular inter- 
est are those that describe the linear response of the "hydrodynamic fields" (mass, energy, 
and momentum densities) to small spatial perturbations of the homogeneous equilibrium 
state. The terminology, hydrodynamic fields, is due to the fact that these are the variables 
expected to obey the phenomenological hydrodynamic equations on large space and time 
scales. The response functions provide the means to study such hydrodynamic excitations 
starting from their fundamental basis in non-equilibrium statistical mechanics. For exam- 
ple, they provide the formally exact Helfand [5] and Green-Kubo [4, 6] representations for 
transport coefficients. More generally, the response functions describe the broader range of 
excitations on shorter space and time scales as well. One of the most instructive theoretical 
approaches to their evaluation has been kinetic theory, with the greatest progress made for 
the idealized fluid of hard spheres [4, 7]. 

Recently, this linear response approach has been extended to granular fluids [8, 9]. The 
objective here is to demonstrate the application of kinetic theory methods for the evalua- 
tion of the granular response functions. Only the case of smooth, inelastic hard spheres is 
considered both for simplicity and to parallel closely the corresponding developments for 
normal fluids. Such an idealized model still captures the most important features of many 
granular fluids [10]. The usual notion of kinetic theory is a nonlinear equation for the prob- 
ability density in single particle phase space. An advantage of the linear response functions 
is that their kinetic equation is inherently linear. For practical purposes, a Markovian ap- 
proximation to this linear kinetic equation is described, based on the neglect of dynamical 
correlations. The approximations leading to this equation, and the differences between its 
implications for normal and granular fluids, are discussed. The nature of the approximation 
does not a priori assume weak dissipation, low density, or large length and time scales. In 
the elastic limit, it becomes the linear Enskog kinetic equation for the response functions of 
a normal equilibrium fluid. A granular Enskog limit is described here as well, by the further 
neglect of all velocity correlations in the Markov approximation. 

A related set of time correlation functions involving the fluxes of the hydrodynamic fields, 
instead of the fields themselves, determine the Helfand and Green-Kubo representations 
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for the transport coefficients of the Navier-Stokes order hydro dynamic equations [8, 9]. 
These can be evaluated by the same kinetic theory developed for the response functions. 
This is described in detail for the shear viscosity in the Markov approximation, indicating 
how the time dependence of these flux correlation functions can be determined, as well as 
the associated transport coefficient. The remaining Navier-Stokes transport coefficients are 
evaluated in the Appendix D. It is confirmed that, in the linear Enskog approximation, the 
results agree in detail with those obtained from the Chapman-Enskog method to solve the 
nonlinear Enskog equation for the distribution function [11]. 

A similar program has already been carried out in the limited context of dynamics for 
an impurity particle in a granular fluid [12, 13]. In that case, the only hydrodynamic field 
is the impurity particle probability density, and the response function is its autocorrelation 
function. The hydrodynamic equation is a diffusion equation, and the Green-Kubo expres- 
sion for the diffusion coefficient is given by the time integral of the velocity autocorrelation 
function. This was evaluated by kinetic theory in the Enskog approximation, and both the 
correlation function and the transport coefficient were compared with molecular dynam- 
ics simulation data over a wide range of densities and degrees of inelasticity. The results 
provide an instructive characterization of the domain of validity for the Enskog (Enskog- 
Lorentz in this case) kinetic equation, and expose important differences between impurity 
dynamics in normal and granular fluids. The presentation here constitutes an extension of 
that theoretical analysis to the full range of multi-particle mass, energy, and momentum 
transport. 

The origin of a kinetic theory is an exact hierarchy (the BBGKY hierarchy [4, 7]) of 
equations for the reduced few particle representations for a property of interest. A kinetic 
equation is comprised of the first hierarchy equation together with a "closure", expressing 
the solution to the second hierarchy equation as a functional of that for the first. This leads 
to a closed, deterministic equation for the latter which is the kinetic equation. Practical 
methods have been developed for normal fluids based on inversion of cluster expansions and 
partial resummations, as well as more phenomenological estimates. This is the point at 
which kinetic theory confronts the difficult many body problem, and one objective of the 
current work is to motivate a corresponding attention to such details for the granular fluid. 
Only in this way can the qualitative speculations about differences between normal and 
granular fluid be made more precise. An example of such uncertainties is the role of velocity 
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correlations in the construction of the closure. It is well known that velocity correlations 
generated by collective many particle collisions (e.g., ring collisions) are responsible for 
the dominant density dependence of transport coefficients at very high densities, but they 
are relatively unimportant at low to moderate densities where simple spatial correlations 
(e.g., excluded volume effects) are dominant. The latter are incorporated in the Enskog 
approximation for accurate corrections to the predictions of the low density Boltzmann 
equation. Granular fluids introduce a complication to this separation of dynamical velocity 
correlations and static structural correlations, according to the density considered. In these 
systems, there are inherent static velocity correlations, not directly associated with many 
particle dynamics, that are present even at low to moderate densities. The kinetic theory 
for response functions provides an appropriate setting for the study of the quantitative 
importance of these correlations on properties of interest. 

The response functions in the Markovian approximation are expressed in terms of the lin- 
ear generator for dynamics. Questions about the existence and dominance of hydrodynamics 
can be made precise at this point, by asking if the hydrodynamic modes (eigenvalues) appear 
in the spectrum of this operator and if they are the slowest modes. The first part, existence 
of hydrodynamic modes, can be demonstrated at long wavelengths if expected conditions of 
analyticity are satisfied (see Sec. VI below). The second issue of dominance at long times 
can be addressed practically using simplified kinetic models for this generator. Such results 
support the primary assumptions in references [8] and [9], for a derivation of formally exact 
expressions for the hydrodynamic transport coefficients in terms of time correlation func- 
tions. However, the linear kinetic equation applies as well to short space and time scales, 
that are important for non-hydrodynamic response. 

The main points of this analysis are summarized in the last section. 

II. LINEAR RESPONSE FUNCTIONS 

An idealized granular fluid of N smooth, inelastic hard spheres (d = 3) or disks (d = 2) 
of mass m and diameter a is considered. The inelasticity is characterized by a constant 
coefficient of normal restitution a. The properties of interest are the average number density 
n(r,t), the granular temperature T(r,t), and the local flow velocity U(r,t), collectively 
denoted by y(r,t) = {yp{r,t)}. The response to be studied here is the variation of these 
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fields at time t due to a variation in their initial values at time t — 0. The initial conditions 
are spatial variations of a homogeneous reference state, yp(r, 0) = yp,h{fy + 0). In 

the most general case, the response 6yp[r,t\y (0)] depends nonlinearly on the perturbations 
8y(r,0), but if these are small, the linear order gives the dominant response, 

5yp[r,t\y(0)] = J2 [ dr' C^(r - r',t) 6y^r',0), (1) 

7 J 

with the linear response functions defined as 

~5yp [r,t\y 



Cfr (r-r',t) 



(2) 

5j/(0)=0 



fy 7 (r',0) 

Since the reference state is homogeneous, the linear response functions depend on r and r' 
only through their difference. 

The difference between granular and normal fluids occurs already at the level of the 
homogeneous reference state. For normal fluids, this is the equilibrium stationary state and 
all the time dependence of the response functions is due to the spatial perturbations. For 
granular fluids, the homogeneous reference state is inherently time- dependent, even without 
perturbation, due to the "cooling" of inelastic collisions. As a consequence, the temperature 
T h of the homogeneous granular fluid decreases in time according with a cooling law 

^ = -Co[T fc (f)]T fc (*), (3) 

where the cooling rate ( (T h ) is a characteristic function of the homogeneous cooling ref- 
erence state. Thus the relevant response at time t is measured relative to the reference 
homogeneous state at the same time rather than to the initial state. Then, dimensionless 
fields 5y*p are introduced by 

r, *i f fyp ] {Sn 5T 5U } 

{Sy ' } ^\f^)rU:m-W)S' (4) 

where the definition of the reference fields yp t h{t) follows from the second identity and 



v (t) = 



2T h (t) 



n l/2 



(5) 



m 

is a thermal velocity. The dependence of the linear response functions on r — r' suggests 
the utility of a Fourier representation of Eq. (1) that is expressed in the form 
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(k*,s) = J2c;,(k*,s)5y;(k*,0), (6) 



where a tilde over a function denotes its dimensionless Fourier transform, defined by 

f(k*) = r d J dre ik *- r / e f(r) 

= J dr*e* k *- r *f(r), (7) 
with r* = r/£, and dimensionless response functions have been identified as 

c; y (k*, s ) = ^i(0^(fc*,0y 7 ,fc(°) 

= T^ h it) J dre ik *^C p ,(r,t)y lth (0). (8) 

Here k* = kt is a dimensionless wavevector, with I being a characteristic length of the 
system. Moreover, a dimensionless time scale 

(9) 
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has been introduced. This time s has the interpretation of an average collision number per 
particle up to time t, if £ is chosen to be the mean free path of the particles. The remainder 
of this presentation focuses on the dimensionless linear response functions C*^. 

For small k* and large s, the functions 5y^(k*,s) are expected to obey the linearized 
hydrodynamic Navier-Stokes equations, and the response functions correspond to the Green 
functions for the solution to the initial value problem associated with those equations [9]. 
That description is only phenomenological, since it is parameterized by the unknown trans- 
port coefficients. A more complete and exact description for all k* and s is provided by 
non-equilibrium statistical mechanics, as described in references [8] and [9]. Briefly, the 
construction is as follows. 

The hydrodynamic fields yp(r,t) are defined from averages of the microscopic number 
density, energy density, and momentum density over the phase space density p(T,t), repre- 
senting the probability that the positions q r and velocities v r of the particles have specified 
values denoted by T = {xi, . . . ,x N }. The latter is a point in the 2Nd dimensional phase 
space with the notation x r = {q r , v r }. For any specified initial state, p(T, 0), this probability 
density evolves in time according to the Liouville equation 

| + r) P (r,t) = o, (io) 



with 

AT 

d 



N a , N N 

a 1 
dq~ r ~~ 2 

r=l ^ r r=l s+r 



L = J2 v r-^--lJ2J2 T ( x r^s)- (11) 



The time independent operator L is the generator for the hard sphere dynamics, where the 
singular binary collisions are described by 

T(x r ,x s ) = 5(q rs - a)\q rs ■ g rs \ [0 (q rs ■ g rs ) a^b^ 1 - (-q rs ■ g rs )] . (12) 

In this expression, q rs = q r — q s , g rs = v r — v s , is the Heaviside step function, q rs = q rs /g rs , 
and 6" 1 is the substitution operator that replaces the velocities v r , v s by their "precollisional" 
values v",v", 

b; s 1 F(v r ,v s ) = F(v , ;,v':), (13) 

v" = v r - ^ a (q rs ■ g rs ) q rs , (14) 
la 

1 + a 

v*" = v s + — — (q rs ■ g rs ) q rs . (15) 
la 

For an isolated system, instead of the equilibrium state for a normal fluid, there is a 
spacial solution to the Liouville equation of the form 

PhF,t) = [£v (t)}- Nd p* h ({q rs /e,v r /v (t)}) = [£v (t)]- Nd pl(r). (16) 

The dimensionless phase point T* = {x*, . . . ,x* N } is now expressed in terms of the scaled 
positions and velocities x* = {q*.,v*} = {q r /£, v r /v (t)}. This solution depends on the po- 
sitions only through the relative variables q rs and, therefore, it has translational invariance, 
representing a homogeneous state. All of the time dependence occurs through the thermal 
velocity defined by Eqs. (5) and (3), so the Liouville equation becomes for pjj(r*) 

Tpun = o. (l?) 

The operator C is the sum of the original generator for trajectories, now in the dimensionless 
variables, plus a scaling operator representing the time dependenc of vo(t) using Eq. (3), 

t* N a 

2*=£*+!z^-< («) 



2 ^dv; 

r=l ' 



where Co is the dimensionless cooling rate 



G = -777 Co PUOl ■ (19) 

The solution to Eq. (17) will be referred to as the homogeneous cooling state (HCS). 
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For more general solutions to the Liouville equation, it is useful to introduce the same 
scaled velocities to account for this inherent cooling. Then, the Liouville equation in dimen- 
sionless form becomes 

(^+r) P *(r*, s ) = o. (20) 

In this form, it is seen that the HCS is a stationary solution to the dimensionless Liouville 
equation. This is an important result for the representation of granular response functions. It 
shows that, in the appropriate dimensionless form, the reference state is again stationary, just 
as for equilibrium fluids. However, the introduction of this stationary representation comes 
at the price of changing the generator for the dynamics from L* to C* . For the purposes of 
the discussion here, it is assumed that all properties of this homogeneous reference state are 
known. Further comments on the HCS are given in the Appendix A. 

The deviations of the relative hydrodynamic fields form their values in the HCS, 
5y*p(k*, s), are the averages of associated phase functions a*p(T*, s), 

5y;(k*, s ) = J dr*a;(r*;k*)\p*(r*, s )- P * h (r*)}, (21) 

with 

1 N 

a;(r*;fc*) = -^e^^a;«), (22) 

where n* h = n h £ d is the number density of the HCS in the reduced units and the single 
particle functions a*p(v* r ) are defined by 

f 2v* 2 1 
= |l, — (23) 

For later convenience, the components of the flow field and the velocity of the particles have 
been chosen to be a longitudinal component along k*, and d — 1 transverse components, 
so that = k • v* and v* u = e, • v*, with jfc = k*/k*, e*; i — 1, . . . , d — 1 j forming a set 
of d pairwise perpendicular unit vectors. To formulate the linear response problem, the 
initial state p*(r*,0) for the solution to the Lioville equation in (20) is chosen to be close 
to the HCS, in the sense that it is a functional of the initial fields 5y%(r*,0) and becomes 
the HCS for Sy^ — > 0. More specifically, the system is viewed as partitioned into small cells 
such that the distribution function is the HCS in each cell, but with different values for the 
hydrodynamic fields. This is the analogue of the local equilibrium distribution for normal 
fluids, and will be referred to as the local HCS, Pi h [T*\5y*]. Its construction is discussed in 
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Appendix A. For small 5y*p, the solution to the Liouville equation for the initial local HCS 
expanded to first order is 



sC* _* /-p* 

KJJ 



p*(r,s)-p* h (r*) = e~ sL P *(r*,o) 

d+2 

^;(r*,0). (24) 

<5j/*=0 



7=1 



<fy 7 (r*) 

Substitution into Eq. (21) allows identification of the response functions C^(r*, s) and their 
equivalent Fourier representation (using translational invariance of the functional derivative 
at 8y* = 0), 

C* M (k, s) = J dT*a; (r*; k*) (r*; 0*, s) , (25) 



with 



and 



^;(r*;r*,s) = e" s£ ^(P;r*) (26) 



^ (r*;r*) = 



(27) 

<5j/=0 



5y* (r*) 

These results are closely analogous to those for a normal fluid, where p| ft [r*|5|/*] is an initial 
local equilibrium ensemble and p^(T*) is the reference equilibrium Gibbs ensemble, p*(T*). 
For example, if the grand canonical ensemble were used, the ^*'s would become linear 
combinations of the set {a*} times the equilibrium ensemble, and the response functions 
would be equilibrium time correlation functions for the local conserved quantities. 

III. KINETIC THEORY 

All of the following analysis is carried out in terms of the dimensionless variables, so the 
asterisk will be left implicit for simplicity. The response functions of Eq. (25) are expressed 
in terms of the full iV-particle phase space. A reduced description in terms of the single 
particle phase space is possible because the functions ap (r, k) defined in Eq. (22) are sums 
of single particle functions. Consequently, integrating over the positions and momenta for 
all except one particle leads to the exact alternative representation: 

dfa(k, s) = — / duia^vi)^ 1 ^!,*;^), (28) 



where 



= J d qi e lk ^\x l)S ). (29) 



The function i(;^\xi,s) is the first element of a hierarchy of iV functions 
ip^ n \xi,x 2 , • • • , x m , s), m — 1, . . . , N, defined through 

~Sf( m h 



ipL ^ {x\ , . . . , x m , s) 



^ (xi , . . . , x m , s) — 



(Xl, 



, Xr 



<fy 7 (O,0) 



Sy=0 



dx m+1 ..dx N p(T, s), 



(30) 



(31) 



(N — m)\ 

1 < m < N . The functions f^ m \xi, . . . , x m , s) are the reduced distribution functions associ- 
ated with the solution to the Liouville equation (20). They obey the corresponding BBGKY 
hierarchy of equations [4] . The first equation of this hierarchy is 



d_ 

ds 1 



Vw)+l^-[t>i/ (i W)] 



dx 2 T{x l ,x 2 )f< 2 \x l ,x 2 ,s). (32) 



dqi J 2 dvi 

Then, it follows directly from the definition (30) that the function ^\xi,s) obeys the 
analogous equation 



dg . dq , i>W(x 1 ,s) + [v 1 ^ 1 \x 1 ,s)} = I dx 2 T(x 1 ,x 2 )^ 2 \x 1 ,x 2 ,s). (33) 



d d 
+ v 1 ■ 



The representation given by Eq. (28) is very appealing, since the iV-particle problem has 
been expressed without approximation in terms of the effective dynamics in the single par- 
ticle phase space. The fundamental difficulty, however, is that the first hierarchy equation 
(33) does not determine this effective dynamics without specifying ip\ '(xi,x 2 , s). An equa- 
tion similar to (33) can be written for ip^\xi,x 2 , s), the second BBGKY hierarchy equation, 
but it in turn requires specification of ip\ (xi, X2, x^, s). In this way, a coupling to the full 
N particle problem recurs. This coupling is broken if i[)j\xi,x 2 ,s) can be specified as an 
explicit functional of ipy\xi,s). It is argued in Appendix B that this functional, when it 
exists, is independent of Sy 7 (r, 0), linear in xjjj\x, s), and has the general form 

^\xi,x 2 , s) = J dxK(xi,x 2 ,s;x)ip^\x, s). (34) 



The kernel defining the functional is 

K(xi,x 2 , s;x) 



5f ( - 2) {x l ,x 2 ,; 
5fW(x,s) 



(35) 



8y=0 



Once K(xi,x 2 ,s;x) is known, Eq. (33) becomes a closed, deterministic, linear kinetic equa- 
tion for ip^ in its most general form, 



d 

OS 



_d_ 
dqi 



+ M(s) 



0(x 1 ,s) + 



2 dvi 



[v^ixus)] =0, 



(36) 



10 



with the formal "collision operator" M given by 

(37) 

This notion of a kinetic equation makes no a priori assumptions regarding the density or 
degree of dissipation, as is sometimes assumed. In fact, it is formally exact so any restrictions 
arise only when specific approximations are introduced to construct the functional. 

The utility of the kinetic theory representation in any specific application depends on 
an accurate construction of the kernel K(xi,x 2 ,s;x) defined in Eq. (35). This is the point 
at which the full many-body problem must be confronted. The idea originated with Bo- 
goliubov [14], who conceived that all f^ m \xi, . . . ,x m , s) with m > 1 take a simpler form, 
/( m ) [xi, ..,x m \f^(s)], after a brief "synchronization" time, where all time dependence occurs 
only through f^(s). For normal fluids at low density, construction of this functional can 
be accomplished by formal density expansions, leading to a sequence of contributions from 
clusters of particles of increasing size [15]. At lowest order, the closure is that associated with 
the Boltzmann equation, while at next order three particle scattering is described. Even for 
normal fluids, these sophisticated cluster expansions have limited direct use due to many- 
particle recollisions events that contribute secular terms in the formal expansions, violating 
the notion of a short synchronization time [16]. These recollisions ("rings") in turn signal 
non-analytic density dependence and slow algebraic decay of correlations in time [17]. The 
experience gained in such studies over the past forty years has provided important insight for 
the construction of more phenomenological closures, such as the Enskog approximation for 
moderately dense gases and mode coupling models for very dense and metastable (glassy) 
fluids [18]. The development of a kinetic theory for granular fluids provides an opportunity 
to revisit many of these issues in an even more challenging context [19]. The analysis of the 
next sections illustrates this for the simplest "mean field" approximation to describe density 
effects beyond the Boltzmann limit. 

IV. THE MARKOVIAN APPROXIMATION 

The approximation described in this section is based on the assumption that the form 
of the correlations and their effect on collisional properties are essentially the same at all 
times. If so, the kernel K(x±, x 2 , s; x) that determines these properties at the two-particle 
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M{x u s)^\ 



x\, s) = — / dx 



dx 2 T(x 1 ,x 2 )K(xi,x 2 , s; x) 



level, can be represented approximately by its form at the initial time s = 0, 



K(x 1 ,x 2 ,s;x) ~ K(xi,x 2 ,0;x). 



(38) 



The collision operator term in Eq. (37) then becomes 

S fth [xi,x 2 \Sy] 



M( Xl )^ (x 1 ,s) 



dx / dx 2 T(xi, x 2 ) 



(i) 
eh 



x 



(i), 



X,S) 



(39) 



8y=0 



where it has been used that the initial conditions f^(xi, .., x m , 0) are the corresponding 
reduced distribution functions associated with the local HCS, p lh [r|5y], given in Eq. (A8), 



fl m Kx 1 ,..,x m ,0)=f™[x 1 ,..,x m \6y]) 



dx m+1 ..dx N p lh [r|5y] 



(40) 



(N — m)\ 

Since the approximate operator M(x±) is time- independent, the entire generator for the 
dynamics of ijffl(x\,s) in Eq. (36) also is independent of time. The equation is written in 
compact form as 



9 . / 



with the generator A(xi) given by 



A(x 1 )X(x 1 ) = 



_d_ 



+ M(xi) 



(41) 



(42) 



for arbitrary X(xi). This is a necessary condition for a Markovian description and, conse- 
quently, Eq. (38) will be referred to as the Markovian approximation. The resulting kinetic 
theory is exact at asymptotically short times, and the nature of the approximation makes no 
explicit limitation on the density or the degree of dissipation. Of course, two particle corre- 
lations that develop over time are neglected. For a normal fluid, this Markov approximation 
leads to the Enskog approximation, where only time independent two-particle correlations 
are taken into account. The neglected time-dependent correlations are found in that case 
to be important only at high densities, and the Enskog approximation provides relevant 
corrections to the Boltzmann results up to moderate densities. It is reasonable to expect a 
similar context for granular fluids, although conditioned by the additional parameter space 
of the coefficient of restitution a. 

The construction of K(xi, x 2 , 0; x) now involves only analysis of the initial reduced dis- 
tribution functions associated with the local HCS, and the corresponding collision operator 
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M is entirely characterized by properties of this state. It is useful to make this more explicit 
in terms of the pair correlation function gf^ [x± , x 2 \ 8y] for the HCS, defined by 



fih [xi , x 2 1 6y] = f$ [x x 1 6y] f$ [x 2 \5y] g$ [xi ,x 2 \Sy}. 

Then Eq. (35) in the Markovian approximation yields 

K(x 1 ,x 2 ,0;x) = g^\x 1 ,x 2 ) f ( h 1) (v 1 )5(x-x 2 ) + fl l 1) (v 2 )5(x-x 1 ) 

6g$[x 1 ,x 2 \6y] 



(43) 



+/i 1) («i)/n«2) 



8fS?[*\8v] 



(44) 



5y=0 



The first term on the right side is given explicitly in terms of the HCS correlations, while 
the second term requires further analysis of the dependence of local HCS correlations on 
fjj^. For a normal fluid, these correlations are independent of the velocities and depend on 
the local equilibrium distribution f^} only through the local density, 

5g$ [xi,x 2 \5y] _^ I 

n h 



<5j/=0 



x ( 2 )r 

09te [91.92 



Sn(q) 



(45) 



8y=0 



( 2 ) 

The functional form of the local equilibrium pair correlation functional g\ e [qi, q 2 \Sy\ is 
known, so explicit construction of the collision operator M is possible in this case. The 
result is the linear revised Enskog kinetic equation [20]. For granular fluids, g eh [xi, x 2 \5y] 
has a more general functional functional dependence on fj^ through its explicit functional 
dependence on 5y, since 



(9) = — / dxi o 7 (xi) S(q- qi) - /^(ui)} , 



(46) 



where the a 7 (xi)'s are the single particle functions given in Eq. (23). Thus 



6g$[xi,x 2 \6y] 
^ Sfi 1} [x\5y] 



^(9) 



5 y =o a L ] Sy=0 

This provides the practical route for constructing the Markovian kinetic theory for granular 
fluids. The collision operator becomes 

M(xi)i#>(si,a) = - j dx 2 T(x 1 ,x 2 )g^(x 1 ,x 2 ) [/^fa^a*, s) + ff\v 2 )^\ Xl , s)] 
~E / dq2Cx(v h q 12 )— [ dv 2 a x (v 2 )^\x 2 ,s), (48) 
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with 



cx(v u q 12 )= J dxT(x 1 ,x)fi 1 \v 1 )fi 1 \, 



5gf^[xi,x\5y\ 



(49) 



5y=0 



Further discussion and simplification of c\(vi t qi 2 ) is given in Appendix C. 

The response functions of Eq. (28) are given in terms of ip^\vi, —k, s), that is propor- 
tional to the Fourier transform of ip^\xi, s) as indicated in Eq. (29). An equation for the 
latter can be easily derived from Eq. (36), 



ij>W(v 1 ,-k,s) = 0, 



(50) 



where the generator for the dynamics is the Fourier transform of A defined in Eq. (42), 



A(fc) = —ik ■ Vi + M (k) + — ( d + v ■ 



(51) 



with M(vi, k) defined by 

MivuQXM = - J dx 2 T(x 1 ,x 2 )g < ^\x 1 ,x 2 ) [e*^tf ) ^ 1 )X( V2 ) + / J [ 1) (u 2 )X(« 1 ) 

- Vc A («ifc)— [ dv 2 a x (v 2 )X(v 2 ), (52) 

for arbitrary X(v±), c\(vik) being the Fourier transform of c\(vi, q\ 2 ). 
The response functions in Eq. (28) become 

1 



Cfa (fe, s) 



(53) 



(54) 



Sy=0 



Equation (53) is the primary practical result of our analysis here. It provides a realistic 
kinetic theory description of the most fundamental time dependent fluctuations in a granular 
fluid: those induced by perturbations of the hydrodynamic fields. To appreciate the scope 
and generality of this result, note that the Markovian approximation is exact at short times 
for all k, densities, and degrees of restitution. In this short time limit, Eq. (53) yields 



IimC7^(fc, a ) = Jm( e -^*))^, 



'07 ■ 



with 



Nfrik) = — I dv 1 a f3 (v 1 )A(k)^ 1 \v 1 ,k). 



(55) 



(56) 
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At longer times, the Markovian approximation is expected to continue to provide a good 
approximation across this parameter space, since the derivation does not explicitly require 
any limitations on k, densities, or degrees of restitution. This expectation is borne out in 
the elastic limit, where the linear Enskog theory for equilibrium time correlation functions 
is recovered, as discussed in the next section. In that case, comparisons with molecular 
dynamics simulations (wavevector dependent transport [21]) and neutron scattering experi- 
ments [22], confirm the accuracy and practical utility of this kinetic theory over a wide range 
of wavevectors and densities. A determination of the corresponding domain of accuracy for 
granular fluids, awaits similar comparisons with simulation and experiments, but there is no 
simple reason to expect qualitative rather than quantitative differences from normal fluids. 

In Sec. VI, it is shown that the operator A(k) contains the hydrodynamic modes in its 
spectrum for small k. Thus the response functions in Eq. (53) provide a means to study the 
transition from short time dynamics to a presumed dominant hydrodynamics at long times 
for granular fluids, in a manner similar to that done for normal fluids [4] . In addition, since 
it is valid for all k, the nature of hydrodynamics beyond the Navier-Stokes approximation 
can be studied. 

An important application of general linear response methods, is the derivation of Helfand 
and Green-Kubo expressions for the transport coefficients [8, 9]. These are formally exact 
results given in terms of time correlation functions. The above analysis for the response 
functions applies to these as well, and their dynamics in the Markovian approximation is 
generated by the same operator A(k). The evaluation of the Helfand and Green-Kubo 
expressions for the shear viscosity is illustrated in Sec. VII. The results provide a general- 
ization of those from Enskog kinetic theory [11], to include pair velocity correlations. When 
such correlations are neglected, the results of ref. [11] are recovered in detail. This is verified 
for the other Navier-Stokes transport coefficients as well in Appendix E. 

V. GRANULAR ENSKOG APPROXIMATION 

The Markovian approximation discussed in the previous section, requires specification of 

(2) (2) 

g^ h (xi,X2) and C\(xi,qi2), or equivalently 5g ih [x±, x\5y]/5y\(q2) for Sy = 0. While these 
are well defined in terms of the local HCS distribution, little is know about their detailed 
forms as yet, except in the elastic limit where they are accurately determined from liquid 
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state theory of the pair correlation function. The important simplification in that case is 
the absence of velocity correlations. It is plausible to assume that such correlations remain 
weak for the granular fluid as well, and to make the approximation 



g^[x 1: x 2 \Sy] ~ g { £[q l ,q 2 \8y] = g^[q 1: q 2 \Sn}. 



(57) 



The last equality recognizes that the neglect of velocity correlations leads to a functional 
that is independent of 5T and 5U (for hard spheres or disks) and hence is a functional only 
of the density. Furthermore, g^[qi,q 2 \0] = <7^(?i2) as a consequence of fluid symmetry. 
Then, Eq. (49) reduces to 

8geh{qi,q\6n) 



(v 1 ,q 12 )=6 xlC ( E \v 1 „ qi2 )=5 xl J dxTix^fPWfPiv) 



Sn(q 2 ) 



8n=0 



(58) 



The collision operator (52) now simplifies to 
M E X{ Vl ) = -gf\a) I dx 2 T( Xl ,x 2 ) 



ik qi2 



-—^ E \ Vl ,k) [ dv 2 X(v 2 ). 
n h J 



(59) 



It remains to give the explicit density dependence for gf^[ri,r\5n] [23]. As a practical 
matter, it can be chosen to be the pair distribution for a nonuniform normal fluid for 
which a well-developed theory exists. In that case, g h (|r*i — r\) = g eq (\ri — r\), the radial 
distribution function for a uniform hard sphere fluid. Also, for this choice the functional 
derivative appearing in the expression of c^- E \vi, k) can be evaluated for the first few terms 
of a k expansion, as is required for evaluation of transport coefficients. The context of such 
a choice would be that the static spatial correlations of a hard sphere system are due to 
excluded volume effects, and these can be captured using the pair correlation function of a 
fluid of elastic hard spheres. The generator A(fc) obtained in this approximation from Eq. 
(59) gives the granular Enskog kinetic theory, and is the linearized version of the one studied 
in ref. [11]. 



VI. HYDRODYNAMIC MODES 



An important feature of the response functions considered here is their relationship to 
hydrodynamic response. At small k and large s, these response functions should correspond 
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to those from the phenomenological hydrodynamic equations. For example, it is this rela- 
tionship that allows the derivation of Helfand and Green-Kubo expressions for the transport 
coefficients. Any acceptable approximate kinetic theory for the response functions should 
preserve this relationship to hydrodynamics. More specifically, the hydrodynamic excita- 
tions should appear in the spectrum of the linear operator A(fc), 

k(k,v)^\k,v) = \ {/3) {k)(f) i/3 \k,v), (3 = l,...,d + 2, (60) 

where the set (A^(fe)} are the eigenvalues of the d + 2 linearized hydrodynamic equations. 
The above eigenfunctions and eigenvalues are determined in the limit k = in Appendix D 
with the results 

{A^(0)} = {0,|,-|}, (61) 

{0 {/3) KO)} = |^ 1) KO)-2^^^ & 1 \v,0),& 1 \v,0),$ 1 \v,0)\. (62) 

The eigenvalue —Co/2 is d-fold degenerate, and the associated eigenfunctions are the com- 
ponents of the vector ip^\v, 0) = —df^\v)/dv. 

In the elastic limit, these eigenvalues are all zero, corresponding to the d + 2 conservation 
laws, and the eigenfunctions become Maxwellians times linear combinations of the summa- 
tional invariants (1, v 2 , v). For inelastic collisions, the nonzero eigenvalues describe response 
of the cooling temperature to linear perturbations and growth of a constant velocity per- 
turbation relative to the characteristic cooling thermal velocity. In both cases, these are 
also the eigenvalues of the phenomenological linearized hydrodynamic equations in the long 
wavelength limit. 

With the eigenfunctions and eigenvalues known at k — 0, their values for finite but 
small k can be obtained by perturbation theory. In this way the Navier-Stokes transport 
coefficients can be determined directly from the coefficients up through order k 2 . This direct 
calculation of the spectrum for the generator of a linear kinetic theory has been described in 
detail recently for the granular Boltzmann equation [24], and its extension to the Markovian 
kinetic theory given here is straightforward. Instead, the remainder of this presentation 
addresses the calculation of the transport coefficients from an approximate evaluation of 
their Helfand and Green-Kubo representations that have been obtained from linear response. 
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VII. HELFAND AND GREEN-KUBO EXPRESSIONS 



As noted above, the exact response functions defined in Sec. II must agree with those of 
the linearized phenomenological hydrodynamic equations in the long wavelength and long 
time limit. This relationship allows identification of the parameters of those phenomenolog- 
ical equations in terms of the response functions in this limit. The results of this analysis 
for granular fluids has been given recently, leading to expressions for the transport coeffi- 
cients in terms of certain time correlation functions derived from the response functions [9] . 
These correlation functions can be evaluated approximately by the Markov kinetic theory, to 
obtain explicit results for all transport coefficients appearing in the Navier-Stokes hydrody- 
namic equations. Further, neglecting the velocity correlations in the Markov theory, allows 
the evaluation of these quantities in the granular Enskog theory, reproducing the results 
reported in [11]. In this section, only the shear viscosity is considered as an example, while 
all remaining transport coefficients are analyzed in Appendix E. 

The exact Helfand and Green- Kubo expressions for the shear viscosity are (dimensionless 
units are still assumed) 

r] = lim % (s) = Vl v H (0) + lim f ds' n v G (s') , (63) 

Jo 

respectively, where the symbol lim denotes the hydrodynamic limit of V — > oo, followed by 
s — > oo. The correlation functions in the above equation are defined by 

= - j2 ^_ 2 EE / dTH^e-^^M^in (64) 
i=i j=i 

Ql(s) = ^-ni(s) = -^^EE / dTH tJ (T)e-<^)T rhlJ (T). (65) 

i=i j=i 

Here, Hij(T) is the volume integrated momentum flux, 

N N N 

H ij( T ) = E Vr > iVr >i + E E H i? Xr "> ' ( 66 ) 

r=l r=l Sy^r 



r( 2 ) ( T rl _(l+ a Vr , _ww : . ^ . ,2 

ij -l>s) — ^ 

ij is the traceless tensor 



H ij ( x r, x a ) = 5 (q rs - a) 6 (-q rs • g rs ) (q rs • 9\ m f qrstfrsj, (67) 
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and T V:i j(T) is the associated Green-Kubo conjugate flux, 

r v , lJ (T) = -(c+ C ^jM v , J (T). (69) 

It is seen in Eqs. (64) and (65) that the generator of dynamics is £ + y. This reflects the 
fact that the k = mode A = —Co/2, of Eq. (61), has been subtracted out. 

The time independent contribution Q V H (0) in the Green-Kubo expression, can be evalu- 
ated exactly from the definitions (66) and (68) with the result: 

n * (o) = e £ / K{vi) \ ( gM ^b + ^-^b - 1 5 ** • ^) ti ){vi) > 

i=i j=i J v l >j L >y 1 / 

(70) 

where 

hij(v) — ViVj + J dxi J dx 2 H < lf (x 1 , x 2 ) K(xi, x 2 ,0; x). (71) 

The first term in the above expression of hij gives no contribution to Q V H (0), from fluid 
symmetry. The second term can be recognized as being proportional to the average collision 
frequency, v av , as determined by the loss part of the right hand side of the hard sphere 
BBGKY hierarchy (32) specialized for the HCS, 

^=wM»>- < 72 > 

vav = 2c/- 1 J da J dvt J dv 2 (-a ■ g 12 ) (a ■ g 12 ) tf ] (a, v u v 2 ) . (73) 
A. Evaluation in the Markov Approximation 

A complete evaluation of the correlation function Q V H (s) is possible using the Markov 
kinetic theory. As in Sec. Ill, the correlation function can be given a representation in 
terms of one and two particle functions 

d d 



Wh(*) = - ^ E£ / dx.v^M^^s) 

i=l j=l 
_^ d d p p 

Tz^£Y.j dx ^j dx * H if ^) m vIj ^ s ) > ( 74 ) 



d? + d - 2 

i=l j=l 



with the reduced functions M-^fj (xi, ■■jXm, s) defined by 

M<$ (x u ...,x m ,s)= {N m m)l J dx m+1 ... J dx N e~< X+ ^)M r ,, 3 (T). (75) 
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Similarly to the functions t/> 7 considered in Sec. Ill, the above functions obey a BBGKY 
hierarchy, the first equation of which is 

(l + f + "-4r)^») + Tsr-[^(»,»)] 

= J dx 2 T(x 1 ,x 2 )M ( ^J j (x 1 ,x 2 ,s). (76) 

The Markovian approximation in the present case is the same as that defined by Eqs. (34) 
and (38) 

M { ^ j (x 1 ,x 2 , S ) ~ j dxK(x 1 ,x 2 ,0;x)M { ^ j (x,s). (77) 
Then, Eq. (76) becomes the Markovian kinetic equation 

| + | + a) M« = 0, (78) 

where the linear operator A is the same as defined in Eq. (42). 

In this approximation, the Helfand expression of the shear viscosity of a hard sphere or 
disk granular fluid becomes 

i=l j=l 

4(^ + ^-3^-^r)^W' (79) 

where hij(v) is defined in Eq. (71). 

Next, the Green-Kubo expression for 77 in the Markov approximation can be identified 
from Eqs. (63) and (65), 

1=1 J=l ,J ' y 

(80) 

and 

^ G W = rf2 2 EE / d*i N,-(«i)e-( A+ ^) 7 « («i) • (81) 

i=i j=i 

The reduced conjugate flux 7^ is 

7« (-) = - (A + f ) i + ^ ~ ■ £) /?W (82) 

Comparison of Eqs. (70) and (80) shows that (0) is given exactly in the Markov approx- 
imation. The Green-Kubo representation for the shear viscosity requires the large s limit of 
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the integral over s in Eq. (63). It can be verified that Qq (s) has no invariant part, so that 
this limit is expected to exist. This issue is discussed in some detail in Appendix E. Then 
the Green-Kubo expression for shear viscosity can be written as 

d d p 

1 = Q » (°) + d^ + d-2 E E / d^KMVn M , (83) 
+ i=i j=i J 

where (v) is a solution to the integral equation 

Upon writing the above equation, it has been taken into account that terms involving spatial 
derivatives give a vanishing contribution to the expression of the shear viscosity. This is the 
traditional form in which expressions for transport coefficients are obtained from a Chapman 
Enskog expansion of a normal solution to the kinetic equation governing the dynamics of 
the system. 



B. Evaluation in the Granular Enskog Approximation 

The further neglect of velocity correlations in the collision operator M, leads to the 
granular Enskog approximation, i.e., the results given by Eqs. (79) and (81) apply with only 
the replacement A by A E , with 

A^ x )^ v .^ + M E (x) + ^.v, (85) 

and the operator M E given by, 

M E X( Xl ) = -g 2 h (cr) J dx 2 T( Xl ,x 2 ) [fj l 1 \v 1 )X(x 2 ) + X ( Xl ) (v 2 ) 

— — / dq 2 c 1 (x l! ,q 2 ) [ dv 2 X(x 2 ). (86) 
n h J J 

The function c\ is given in Eq. (58). In the Enskog approximation, the conjugate flux 7^ in 
Eq. (82) becomes 

e . Co \ 1 / 9 , 8 2, a \ 



1, M — (A^f)- 2 ^ + ^ " ~M ■ Wi ) f?M. (87, 

One final simplification occurs for the shear viscosity and some other transport coeffi- 
cients. The mean field term in Eq. (86) vanishes when acting on 7^ (v) and, therefore, 

(A* + |) Hj (v) = (j(x) + |) 7 « (*) , (88) 
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where the operator J{x) has been introduced, 



(89) 



Here, X is the linearized Boltzmann collision operator for inelastic hard spheres or disks, 

JX(xi) = J dx 2 T(x u x 2 ) [fl 1 \v 1 )X(x 2 ) +X(x 1 )fi 1 \v 2 )] . (90) 

Therefore, the correlation function in Eq. (81) and the expression for the shear viscosity in 
Eq. (83) take the final forms 

{l + a)a d n h7 r d / 2 g£\ay 



i 

<f + d - 2 

d d 
i=l j=l 



1 + 



4T(^) 



f iVj exp 



-5 J + 



Co 



7u» 



and 



<(o) + E E / rfv h ^ v 5 (*) > 
^ i=i j=i j 



(91) 



(92) 



respectively. The Enskog approximation for (0) is 

H 1 ; 2(rf 2 + 2rf)r(^) y 

and is a solution to the integral equation 



^ki-^i/^w/r^) (93) 



(1) 



(^ + |)^5(«)=7«(«), 



(94) 



The above result for the shear viscosity agrees in detail with that obtained in ref. [11], 
through a Chapman-Enskog procedure applied to the non-linear granular Enskog equation. 

The results in Eqs. (81) and (91) for the Green-Kubo integrand are new. At the formally 
exact level, the integrand is given by the correlation between the flux and the conjugate flux. 
In detail, the contributions from and the T(i, j) terms of C, appear to yield singularities 
at t — 0, signaling a possible nonanalytic dependence on t. This occurs even in the elastic 
limit, and is a peculiarity of hard particle dynamics. Consequently, previous theoretical and 
simulation studies have avoided this by studying the Helfand forms for transport coefficients. 
The approximate kinetic theory described here gives an explicit analytic estimate for Vt G (s) , 
whose integral yields a good estimate for the transport coefficients. This suggests that Vt G (s) 
may have a dominant analytic part with a relatively small non-analytic correction. 
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VIII. DISCUSSION 



Kinetic theory has been used extensively as a formal tool for approximate evaluation 
of response functions, in the study of hard spheres as a model for normal fluids. The 
objective of this work is to take a first step in the development and application of this 
tool in the analogous field of granular fluids. The advantage of developing kinetic theory 
in the context of linear response functions, lies in the fact that the resulting theories are 
inherently linear, and provide a more tractable setting to explore questions such as aging 
to hydrodynamics and short wavelength behavior of the exact hydrodynamic response. The 
two primary contributions here are: 1) the development of a practical kinetic theory for 
an important class of granular time correlation functions and, 2) the demonstration of its 
utility for the evaluation of Helfand and Green-Kubo expressions for Navier-Stokes order 
transport coefficients. 

The linear kinetic theory is summarized by Eq. (50). Based on corresponding studies for 
the elastic limit of this equation, it is expected to have a wide domain of validity with respect 
to space and time scales, as well as densities. The nature of the approximation, short time 
functional relationship, does not explicitly entail questions of inelasticity so it is expected 
to apply as well for a finite range of inelasticity. It encompasses the granular Boltzmann 
equation, in the low density limit, and the familiar Enskog equation in the elastic limit. 
The focus here has been on hydrodynamic response, but the theory includes hydrodynamics 
beyond the Navier-Stokes approximation, and even describes very short wavelength non- 
hydrodynamic behavior, that can be more important at moderate and high densities for 
granular fluids. Finally, this kinetic theory applies beyond the set of hydrodynamic fields 
considered here. For any observable, z(T; r), that can be written as a sum of single particle 
functions so that 

N 

z(r,k) = Y,e ih ' qr *M, (95) 

r=l 

and for the same initial perturbation as considered in Sec. II, the response functions are 
given in the appropriate units by 

C~(k;s) = — I ' dvz{v)e- sl ^i)^\v,k). (96) 
n h J 

This opens the possibility to study a wide range of experimental probes and also fundamental 
questions such as the relationship between fluctuations and response. 
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The second contribution, evaluation of the formal representations for Navier-Stokes trans- 
port coefficients, begins the process of exploring the utility of such formal representations, 
as well as verifying their consistency with earlier Chapman-Enskog based studies. These are 
long wavelength properties of the kinetic theory, and therefore a more controlled context for 
its tests. For example, the Markovian approximation provides a practical context for the in- 
troduction of velocity correlations associated with the reference homogeneous state [25, 26]. 
Their effect on transport coefficients at strong dissipation is expected to be important but 
has not been quantified to date. The kinetic theory scheme developed here also provides the 
basis for formulating and assessing more complex theories, such as those describing mode 
coupling dynamical correlations, which are expected to dominate at very high densities. Fi- 
nally, the verification of the agreement between the kinetic theory evaluation of the Helfand 
and Green-Kubo representations here and the earlier Chapman-Enskog method provides 
further support for the implicit assumptions of these complementary formal approaches. 

An interesting new result, both for normal and granular fluids, is the expression of the 
Green-Kubo time-correlation function Qg( s )- In the Enskog approximation, the correspond- 
ing function for the shear viscosity is given by Eq. (91). To interpret this result, the time 
dependence may be estimated from a leading order cumulant expansion, 

Clg (s) ~ nf (0) e~ s '\ (97) 

i = Et/^M") (J+jtHAv) (98) 

T Etj I dv hij^j (v) 

The corresponding Helfand correlation functions, can be inferred directly from this, 

< (s) ~ Qf (0) + nf (0) r (1 - e-' r ) , (99) 

and the shear viscosity can be identified as 

77 ~ nf (o) + nf (o) r. (ioo) 

These results expose the qualitative nature of the time dependence in each case. The re- 
sulting shear viscosity in these approximations agrees with that obtained by a leading order 
solution to the integral equation (94) as an expansion in Sonine polynomials [11]. 

In conclusion, it is hoped that this work provides a starting point to explore systematic 
analytic approximations to the hydrodynamic response of a granular fluid, with the same 
attention to detail given in the context of normal fluids. These, together with numerical 
studies of exact results, provide a means to understand transport mechanisms in this system. 
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APPENDIX A: LOCAL HOMOGENEOUS STATE 

The local HCS ensemble chosen as the initial perturbation of the HCS, represents a system 
decomposed into spatial cells, each in a HCS with its own local temperature, density, and 
flow velocity. It is constructed formally as follows. First, the HCS ensemble is determined as 
the solution to the homogeneous, stationary Liouville equation (17) in dimensionless form, 

^'^^({tw^ 1 4 (A1) 

where Uh, Th, and iih (not shown explicitly) are the flow field, temperature, and particle 
number density characterizing the HCS. Next, a conservative external force is added to Eq. 
(12) keeping the same C* operator, 

L r=l 

Here 4>* ext = 4> ext /2T h , with 4> ext {r) being the potential associated with the external force. 
The solution of Eq. (A2) is, therefore, a function of this potential, 

*/ */ f \ Qrs V r -U h (pext(qr) , AT \\ , . Q x 

ft= "liT' ~MiT' ; r ' s = l - ■ ■ • ■ N s) ■ (A3) 

This can be considered as the nonuniform fluid ensemble corresponding to the uniform limit 
p* h , since in general the density will be nonuniform through its functional dependence on 
4>ext (r), 

n = n[r\(t) ext } , (A4) 
<Pext = <Pext [r\n\. (A5) 

The second equality assumes the functional dependence of the density on the external po- 
tential is invertible so that the potential can be expressed as a functional of the density field. 
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d 

dq; 



_d_ 

dv* 



Ph = 0. 



(A2) 



For normal fluids in the equilibrium Gibbs state, density functional theory assures that this 
is the case. In particular, for any chosen density field there is a unique external potential 
creating that field from the uniform state. It will be assumed that these properties hold 
as well here for the granular fluid, so that Eq. (A3) can be expressed in terms of the local 
density instead of the potential, 

*/ „*// / J Qrs v r ~ Uh , 



Ph = Ph [\j> i^f' r,* = l,...,Nj ) . (A6) 

With p* h " known from the solution to Eq. (A2 ), the local HCS is constructed by the replace- 
ments 

v r -U h v r -U h - 5U(q r ) 



v (T h ) v [T h + 5T(q r )} ' 

n{q r ) -> n h + 5n{q r ), (A7) 



to get 



Note that the local HCS is no longer a solution to any Liouville equation, but rather is 
simply a reference ensemble representing an hypothetical HCS with different hydrodynamic 
parameters in each spatial cell of the fluid. Its construction in the way presented above, 
supports that interpretation in the sense that p\ h [r*|5y* = constant] = p* h (T*; y* h + Sy*) 
so that both {p* eh [^*\Sy*]} Sr=0 = p* h (T*;yl) and all functional derivatives of p} h [T*\Sy*] 
become derivatives of p* h (T*;yl) at 5y* = 0. More explicitly, it is 



J dr 1 ...J dr p 



5 p Pih[T\5y] 



Sy a (ri) . ..5y p (r p 



de Ph {T;n h ,T h ,U h ) 
dy a ,h ■ ■ ■ dyp,h 



, (A9) 

u h =o 



y={n h ,T h ,0} 

where T is a point in the phase space associated to the original positions and velocities 
and p(r) the corresponding density. It is instructive to carry out this construction of the 
local ensemble for the case of a normal fluid. Then Eq. (Al) gives the familiar equilibrium 
Gibbs ensemble, and Eq. (A2) gives the same ensemble with the Hamiltonian modified to 
include the external potential. Finally, the construction in Eq. (A8) gives the familiar 
local equilibrium ensemble used in linear response theory for spatial perturbations of the 
equilibrium state. In the grand ensemble the dependence on the local density is implicit 
through a local chemical potential p = p[r\n}. 
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APPENDIX B: TWO PARTICLE FUNCTIONAL 



The reduced distribution functions associated with the solution to the Liouville equation 
(20) are defined as 

iV' f f - 

f (m) (x u ...,x m ,s) = — — — yy j dx m+1 ... j dx N e- sC p eh [T\5y}. (Bl) 

As is done in the main text, the asterisk indicating the use of dimensionless variables is left 
implicit. The /( m )'s are clearly not independent functions. For example, f^(x 1 ,x 2 , s) is 
related to f^(x 1: s) by 

(N-l)fW( Xl ,s) = J dx 2 f^(x u x 2 ,s). (B2) 
This implies that f^(x±,x 2 , s) has the representation 

/< 2 >(*i, x 2 , s) = s)fV(x 2 , s)g^(x 1 ,x 2 , s), (B3) 

where the pair correlation function g( 2 \xi,x 2 , s) has the properties 

g (2 \x 1 ,x 2 ,s) = ^ (2) (x 2 ,x 1? s), (B4) 

J dx 2 f( 1 \x 2 ,s)g ( - 2 \x 1 ,x 2 ,s) = J dx 1 f^(x 1 ,s)g^(x 1 ,x 2 ,s) = N -1. (B5) 

This in turn shows that g^ 2 \xi, x 2 , s) is a functional of f^(x 1: s). Quite generally then, 
f( 2 \x\,x 2 , s) can be considered a functional of f^(xi, s), 

f (2 \x 1 ,x 2 ,s) = fW[x 1 ,x 2 ,s\f (1 Hs)}. (B6) 

However, this functional relationship is not unique. The utility of Eq. (B6) lies in discovering 
a choice that leaves the simplest functional dependence. In the low density limit, where it 
is expected that g^ 2 \x\,x 2 ,s) tends to unity, this is clearly the case, since the functional 
becomes a constant independent of xi,x 2 , and s. More generally, finding an appropriate 
functional form for g( 2 \xi,x 2 , s) requires the detailed analysis of the full many-body prob- 
lem. 
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The corresponding functional for ^\x±,x 2 , s) defined by Eq. (30) can be computed as 



ip^\x 1 ,x 2 , s) 



<fy 7 (0,0) 



Sy=0 



dx 



dx 



tf/< 2 > [x 1 ,x 2 , S \f^(s)} 5fW(x,s) 



5fW(x,s) <fy 7 (0,0) 
*/ (2) [xi,x 2 ,a|/W (s)} 



6fW(x,s) 
= J dxK(x 1 ,x 2 ,s;x)ip! y 1 \x,s), 



Sy=0 

0(x,s) 



Sy=0 



(B7) 



with K(x 1 ,x 2 ,s;x) given by Eq. (35). This entails the additional requirement that 
the functional f^[xi, x 2 , s\-\ is independent of the specific initial fields 5y(r,0), i.e. 

[x 1: x 2 , s|/ <1 - ) (s)] depends on these fields only through f^(x,s). Consequently 
K(xi,x 2 , s;x) is also independent of such initial data. This is expected, since the collision 
operator M(s) constructed from K{x\, x 2 , s; x) by means of Eq. (37) should be universal for 
a wide class of initial conditions. Finally the functional form for ip^\xi,x 2 , s) is seen to be 
linear in ipj\s), while in general f^[x\,x 2 , s\f^ (s)] is a nonlinear functional of (s). 



APPENDIX C: INTERPRETATION OF c x (v,k) 

The contribution to the action of the collision operator M(k) on ipy\vi, —k, s) from the 
term proportional to c\(vi, k) in Eq. (52) is 

- J^CA(wi,fe)— / dva x (v)i;^\v,-k,s) = -^2cx(v u k)Cx y (k;s) , (CI) 
x UhJ X 

where Eq. (28) has been employed. The above contribution depends on only low order 
moments of the dependent variable ip^ in the kinetic equation (50), and in fact only those 
moments are of interest for determining the response functions. In this sense, c\ is a mean 
field operator rather than a true collision operator, since its action does not depend directly 
on differences in ip^ before and after a collision like the first term of (52). Instead, Eq. 
(49) shows that c x reflects an average of collisional effects induced through changes in the 
correlations. To provide some interpretation of this term, consider first the elastic limit. 
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1. Elastic limit 



In this case, g^[xi,x 2 \5y] = gf^[qi^q 2 \8n\^ independent of the velocities, temperature, 
and flow field, i.e. it is a function of the spacial coordinates and a functional of the density. 
Equation (49) becomes 

8gfJ[qi,q\Sn] 



c\(v u q 12 ) = 5 X in h / dxT(x 1 ,x)ip(v 1 ) ip (v) 



5n(q 2 ) 



(C2) 



Sn=0 



where •p (v) is the Maxwellian. The pair correlation function for a nonuniform fluid, 
9e e [9i ' ^2 1 8n] , appears in the stationary first BBGKY hierarchy equation ( 32) in the pres- 
ence of an external potential (f> ext associated with the given density (see Appendix A), 

d ( d , r , A d 



n(qi)¥>(ui) 



= J dx 2 T(x 1 ,x 2 )(p(v 1 )(p(v 2 )n(q 1 )n(q 2 )g^ ) [q 1: q 2 \Sn\ . (C3) 

The functional derivative of this equation with respect to Sn (q 3 ) evaluated at at 5n = 
gives 



d 



dqi 



n h p(v 1 )v 1 - ^—C(q 13 ) = n h I dx 2 T(x 1) x 2 )<p (uj) y? (v 2 ) gf ] (q 12 ) [5 (q 31 ) + 5 (q 32 )} 



+n 2 h J dx 2 T(x 1 ,x 2 )(p(v 1 )(p(v 2 ) 



&9?e [gl 1 g2H 

<Jn(qf 3 ) 



.(C4) 



<5n=0 



where it has been used that g e (q\ 2 ) = ge e [qi, q 2 \0] and c (q± 3 ) is the direct correlation function 
defined by [27] 



n^c (qfi - q 3 ) = 5 (<?i - q 3 ) + 2n* 



5n(q 3 ) 



(C5) 



The first term on the right hand side of Eq. (C4) can be evaluated using the elastic limit of 
the explicit form for T(x\,x 2 ) given in Eq. (13). The result is 



J dx 2 T(x 1 ,x 2 )<. 



P (vi) <p (v 2 ) gf ] (q 12 ) [5 (q 31 ) + 5 (q 32 )} = gf ] (a) (p (vi) S(q 13 - a)q 13 ■ v 1 



(2) 



(C6) 



Finally, therefore, Eq. (C4 gives 

n h / dx 2 T(x 1 ,x 2 )ip(v 1 )ip(v 2 ) 



bgfe [QuQ2\Sn\ 
8n(q 3 ) 



Sn=0 
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= wi • j— [c(q 13 ) - #(a)9(g 13 - a)] , (C7) 

where the delta function in Eq. (C5) has been written in terms of the derivative of ©(913 — 0"). 
Using this into Eq. (C2) gives the desired result, 

d 

ca(«i, 912) = <5ai?w(ui)vi ■ ^— [c(q 12 ) - ^ 2) (°")@(9i2 - o)] . (C8) 

Cqi 

In the elastic case, it is seen that the contribution shown in Eq. (CI) is the same as that 
for an external force whose potential is 0(913) — gj^ (o) (q 13 — a). The direct correlation 
function has a discontinuity at 913 = o, with value c (o) = g e (o), so the subtracted theta 
function contribution assures that this potential is continuous. 



2. Inelastic collisions 



For inelastic collisions, the effects of c\(vik) are more complex and more difficult to 
interpret. However, a significant difference from the elastic case can be seen already for the 
simplest case of k = 0. From Eq. (48), it follows that 

dg { h\xi,x;n h ,T h , U h ) 



ca(«i 0) = / dq 12 c x (x lt q 12 ) = I dxT{x 1 ,x)f^\v 1 )f K h '\v) 



dyx,h 



where use has been made of the identity 

Sg^[x u x 2 



I 



dq 



Syx(q) 



8y=0 



dg { h\xi,x;y) 
dyx 



£4=o 
(C9) 



(CIO) 



{n,T,U}={n h ,T h ,0} 



This k = limit vanishes for elastic collisions, as it can be seen directly from Eq. (C8), but 
is nonzero for inelastic collisions. This can be verified in the Enskog approximation, where 
velocity correlations are neglected and, therefore, 

dg^\a;n) 



c>i,0) = S X1 



dn 



n=n h 



J dxT{x^x)^\ Vl )^\v) 



d In g^\a\ n 
dn 



'XI 



2 8V! 



(Cll) 



n=n h 



The second equality follows from the first hierarchy equation (32) particularized for the 
HCS, 



(o_d_ 

2 dV! 



VifPM = I dxT(x 1 ,x)g):>(x 1 ,x)f^(v 1 )f^(V), 



(C12) 
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when velocity correlations are neglected. Thus, it is seen that c\(v lt ,0) includes changes in 
the correlations of collisional effects associated with cooling. 

More generally, the Enskog approximation for arbitrary k reads 



cf(t>i,fc) = £ A i J dqe ik q c x {x u q) 

= J dx 2 T(x 1 ,x 2 )fl 1 \v 1 )fi 1 \v 2 ) J dqe^ 

APPENDIX D: HYDRODYNAMIC MODES OF A(0) 



Sg ( S[qi,q2\Sy] 
Sn(q) 



. (C13) 



<5n=0 



In this Appendix, some of the details leading to the solution of the eigenvalue problem 
(60) at k = are given. Consider first the functions i(;^\vi,k), defined in Eq. (54), at 
k — 0. By translational invariance 



<fy 7 (0) 



8y=0 



Sy 7 (r) 



and so 



^K,0) = / dq 



S fS?[zi\6y] 



Sy=0 



(Dl) 



8y=0 



dr 



<5y 7 (r) 



5y=0 



dy 1 



(D2) 



V=Vh 



The last equality is a consequence of the construction of the local HCS, assuring that all 
functional derivatives in the homogeneous limit are related with ordinary derivatives of the 
HCS (see Appendix A). Use of the expression of the operator A, Eq. (51), yields 



A(0)#»( t)l) 0) = M(0) 



dfP (vi;y) 

dy~, 



+ 



y=yh 



2 dv! 



Vi 



dfP fay) 

dy 1 



(D3) 



V=Vh 
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Next, Eq. (52) gives 
M(0) 



dfi 1] (v^y) 



y=Vh 



-J dx 2 T(x 1 ,x 2 )g^\x 1 ,x 2 ) ^ 
[ dxT(xi, x)fh\vi)ffr\ 



dg { h\xi,x; 5y) 



y=Vh 



dyx 



y=yh 



x / dv 2 a x (v 2 ) 
d 



dfj? (v 2 ;y) 



dy 



dy 1 



7 

y=yh 



dx 2 T(x u x 2 )g y h >(x u x 2 ;y)f h '{v^y)^ \v 2 ;y) 



y=yh 
(D4) 



where Eq. (C9) has been employed. The right hand side in the above equation can be further 
simplified by means of Eq. (C12), after changing v r into v r — U, 



M(0) 

oy-t 

and Eq. (D3) becomes 

A(0)^)( Vl ,0) 
Realizing that 

and 



y=yh 



d(o(y) 
dy 1 



_d_ 

dy^ 



1 d 



y=yh 



2dv l 



Co(y) d 

2 dvi 



{v 1 -U)f^\v 1 ,y) 



(D5) 



y=Vh 



Co fdU 



2 \dy. 



} f ( u\vi)- (D6) 



~i ' y=Vh 



dv 1 



9fh ] (v,y) 
dy 2 

dti\v,y) 



ld_ 
2dv 



y=yh 



dy 3 



vfP (v) 
dfi 1] (v) 



dv 



(D7) 



(D8) 



y=Vh 



Eq. (D6) is seen to be equivalent to 



A(0)$ 1} KO) 



d( (y) 
dy 1 



Co fdU 



2 \dy y 



</? } K0). (D9) 



y=yh x a ' / y=yh 

Above, the index 3 and vectorial notation has been used to identify the d components 
associated to the velocity field. 

The specific cases following from Eq. (D9) are 

d( 



A(0)ri 1} (t;i,0) =n h 



dn h 



(D10) 
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and 



A(0)^ 1) K,0) = |^ 1) ( Vl ,0), 



A(0)^ 1) K,0) = -|^ 1) («i,0). 



Equations (D10) and (Dll) can be combined to give 
A(0) 



d In n h y T 



0. 



(Dll) 
(D12) 

(D13) 



The above equations (D11)-(D13) have the form of eigenvalue equations 

A(O)0 (/3) (O, Vl ) = \ {/3 \0)^\0,v 1 ), (D14) 
P — 1, . . . , d + 2, with the eigenvalues and eigenfunctions given in Eqs. (61) and (62). 

APPENDIX E: TRANSPORT COEFFICIENTS 



In this appendix, the pressure, the cooling rate, and the Navier-Stokes transport coeffi- 
cients are evaluated in the kinetic theory approximations developed in the main text. These 
quantities are identified from the linearized phenomenological Navier-Stokes equations for 
an isolated granular fluid, that in dimensionless units have the form 

d+2 



d 



-^;(fc*, s ) + ^/c;^;(fc*, s ) = o, 



(El) 



7=1 



where the transport matrix IC* hyd is found to be block diagonal with a "longitudinal" part, 
corresponding to the fields {5n* , 5T* , 5U]_} , given by 

/ n — ^ 







-ik* 



\ 



J 



CoS + C-f - C n ) ^ f + - C T ) k* 2 -i + c") k* 

ipj dlnp h , * ip* h ,* Co* , [ 2(^-1) . *1 u*2 

2 dlnn h K 2 K 2 [ d " ~+~ K \ K 

(E2) 

The "transverse" components, 5U* ±i) decouple from the longitudinal ones, and their trans- 
port matrix reads 

JC hyd *(k) =(-Q + ^k* 2 ] I, (E3) 



[~f + V*k* 2 ) I, 

with I being the unit matrix of dimesnsion d — 1. These expressions include the unspecified 
functions ph(n h ,T h ) and ( h (n h ,T h ), defining the pressure and the cooling rate of the HCS, 
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respectively, as well as the unknown shear viscosity r)(nh,Th), the bulk viscosity K,(nh,Th), 
the thermal conductivity \(rih,Th), and the new coefficient, fi(rih,Th), associated with the 
contribution of the density gradient to the heat transport in a granular fluid. Finally, 
( n (rih, T/j), ( T (nh, T/j), and ( u (rih, T^) are transport coefficients arising from the local cooling 
rate. The dimensionless form of the cooling rate was defined in Eq. (19), and the definitions 
for the remaining dimensionless quantities are 

Ph * V « 



Ph = — ^r, V = j— , « 



n h Th ' mn h lvo mn h lv ' 



ln h v ' /T/j^o ' 

r = c ^ r = !*C: C* T = t^- (E4) 

Formal expressions for these parameters have been obtained in ref. [9], by carrying out 
a linear response analysis. Each of those expressions is evaluated below by means of the 
kinetic theory developed in the text, for comparison to the results obtained by a Chapman- 
Enskog solution to the nonlinear Enskog equation [11]. First, some general considerations 
that apply to all of the transport coefficients will be addressed. Dimensionless units will be 
used in the remaining of this Appendix, and the asterisk will be suppressed for simplicity, 
as done in the main text. 



1. Reduced Form of the Transport Coefficients 



The Helfand form for a generic dimensionless transport coefficient x as obtained in ref. 
[9] is 

X = lim^(s)=hm^ J dT j {V) (1 - V) e~< c - x ) M {Y) , (E5) 

where j (Y) is a flux associated with the densities of the hydrodynamic variables, and A 
denotes one of the eigenvalues in (61). These fluxes have the generic form 



N N N 



(E6) 

r=l r=l s^r 

ji and j 2 being one and two-particle functions of the phase point, respectively. The adjoint 
functions M. (r) are related to the conjugate densities in the hydrodynamic response defined 
in Eq. (27), and read 



M 



(Y) = M j 



drkr 



5pih[Y\y] 
8y(r) 



(E7) 



y=yh 
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M being some constant. Moreover, V is the projection operator defined by 

, d+2 . 

p(r)x(r) = -^^ 7 (r;0) / dr'a 1 (r'-,o)x (V) , (E8) 

7=1 J 

where the phase functions a 7 (r; 0) and *0 7 (r; 0) are the densities defined by Eqs. (22) and 
the Fourier transform of (27), respectively, both evaluated at k — 0. The generator for the 
dynamics is £ — A, where A is one of the hydrodynamic modes at k = identified in Eq. (61). 
Both the projection operator and the additional time dependence of the term containing A 
are necessary to ensure that the long time limit of the correlation function in Eq. (E5) is 
well defined. 

Proceeding as in Sec. VII for the shear viscosity, a reduced expression for these transport 
coefficients in the Markov approximation can be obtained. The generic correlation function 
£l H (s) of Eq- (E5) becomes 

n„(s)^± J dx 1 J(x 1 )[l-V^(x 1 )]e- s ^M^(x 1 ). (E9) 

The direct flux J (xi) in this reduced time correlation function is 

J (x) = ji (x) + J dxi J dx 2 j 2 {xi,x 2 ) K (xi, x 2 , 0;x) , (ElO) 

where K (x± : x 2 , 0; x) is the kernel for the collision operator given by Eq. (44), and is 
the one-particle function in the hierarchy associated with M. 

M {1) { Xl ) = N J dx 2 ... J dx N M{Y). (Ell) 

The generator for the Markov dynamics A is given in (42). Lastly, is the one particle 
analog of V defined in Eq. (E8), 

1 d+2 ~ If 
V { - 1 \x 1 )X{x 1 ) = -Y J ^ ) {^)— / dxa y (v)X(x), (E12) 

V 7=1 Uh J 

where the ipy\vi, 0) are the functions defined in Eq. (D2). 

Equation (E5) can be transformed into the Green-Kubo form to give 

x = Q H (0) + lim / ds' -^Qh (s') = n H (0) + lim f ds' tt G (s 1 ) . (E13) 
Jo & s Jo 

The Markov approximation in this representation follows directly from Eq. (E9), 

n H (0) J d Xl J (xi) [1 - V^\x 1 )} M m (xi) , (E14) 
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n G (s)^yj dx, J ( Xl ) [1 - V m ( Xl )} e^ A ' x h ( Xl ) , (E15) 
with the conjugate flux j(xi) given by 

7(xi) = - (A - \) M (1) ( Xl ) . (E16) 

In this Green-Kubo form, the role of the projection operator can be readily interpreted as 
follows. Taking into account that projects over the subspace spanned by the hydrody- 
namic eigenfunctions of A(0), the property 

[1 - V {1 \ Xl )] e- s(A -^ = [1 - V (1 \ Xl )] e' 3 ^ [1 - V (1 \ Xl )] , (E17) 

is obtained. This shows that the presence of the projection operator in Eq. (E15) ensures 
that the generator of dynamics exp[— s(A — A)] acts on a function that is orthogonal to its 
invariants and, therefore, has a well defined long time limit. 

The above expressions can be specialized to the granular Enskog approximation by replac- 
ing the Markovian kernel K in Eq. (E10) with its Enskog approximation, and by replacing 
the linear operator A by A E , defined in Eq. (85). Further, the reduced conjugate functions 
in (Ell) and Eqs. (E12) in the Enskog approximation become 

{M^( X )} ~ {Mfc ■ qf£\v), y % -^- » ~Mk ■ q^ft\v)} (E18) 

and 

W\v,0)} * [«/«(„)] (E19) 

respectively. Now, f^\v) is the one-particle distribution of the HCS obtained from the 
nonlinear Enskog kinetic equation. In the rest of the appendix, explicit expressions for the 
hydrodynamic parameters are given in the granular Enskog approximation. 



2. Evaluation in the Enskog Approximation 

The Helfand and Green-Kubo expressions for the Navier-Stokes transport coefficients of 
a hard sphere granular fluid are reported in ref . [9] . From those expressions, the N particle 
functions j (T) and M (T), together with the eigenvalue A appearing in Eq. (E5) for each of 
the transport coefficients, can be read off. Then, following the procedure illustrated above 
and, in the context of the shear viscosity in the main text, the results reported below are 
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obtained. Attention is restricted to the transport coefficients that were calculated in ref. [11], 
which excludes the transport coefficients ( T and C n from the local cooling rate. Moreover, 
the parameter £ defining the length scale is chosen such that n h £ 3 = 1 in the following, for 
the sake of simplicity. 



a. The pressure 

The expression for the pressure has been identified in [9] as 
[1 + a) a 



p = 1 + 



2Vd J dx! J dx 2 5(q 12 - cr)Q(-q 12 - g 12 )(qi2- gi2) 2 fh\qi2,v 1 ,v 2 ). (E20) 

This is the second moment of the normal component of the relative velocity averaged over the 
two-particle distribution at contact. In the Enskog approximation, all velocity correlations 
in the two-particle distribution are neglected, and the above expression simplifies to 

ir d / 2 (l +a)a d m . , _ , 

p=1+ mm gi (g) - (E21) 

b. The zeroth order cooling rate (o 

The homogeneous dynamics in Eq. (El) is determined entirely determined by (q, for 
which the expression 

Co = 1 2 y^ J dxx J dx 2 5{q 12 - a)Q(q 12 - g 12 )\q 12 - g 12 \ z f^\q l2 ,v 1 ,v 2 ) (E22) 

was derived in [9] . When the Enskog approximation for the two-particle distribution function 
is used, this simplifies to 

(2) ((T)(1 _ a2) , r 

Co = 2r ^ J j dv ^ J dv z9nr h } M f ( h ] (v 2 ) 

= ^r^J {p ~ 1] I dVl I dV2 ^ {vi) {V2) ■ (E23) 

c. The Euler transport coefficient ( u 

This is a new transport coefficient unique to granular fluids. Its expression, as obtained 
in ref. [9], reads 

( u = lim Q% (s) = (0) + lim I ds' 0% (s') , (E24) 

Jo 
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where 



and 



Q% ( s ) = V- 1 j dYW (r) (1 - V) e'< c+ ^)M c u(Y) (E25) 



ng(s) = V- 1 J dTW(T)(l-V)e-< c+ ^)T (; u(T). (E26) 

In the above expressions, W (T) is the source term in the microscopic energy balance equation 
that characterizes the dissipation due to the inelastic collisions, 

^ 2 N N 

W{T) = —^^^5{q rs - a)Q{-q rs ■ g rs )\q rs ■ g rs \\ (E27) 

r=l s^r 

and the conjugate density and flux are 

N „ 

M c v(r) = -J2<lr--^- PhF), (E28) 

r r 

T c u(r) = -(z+^jM c u{r). (E29) 

The Helfand form for this transport coefficient in the Enskog approximation can be obtained 
by application of Eq. (E9) with the result 

^ (s) ~ V' 1 J dxJ c u(x) [1 -V {1 \x)} e-< AE+ ^)M^(x), (E30) 



where 



(E31) 
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= - Ql • —^(^,0), (E32) 

and A E is the collision operator given in Eq. (85). The Green-Kubo form of this transport 
coefficient is determined by & H (0) and Qq (s). Direct evaluation of the former gives [9] 

n£(0) = -^(l-a)(p-l). (E33) 

Here p is the pressure in the Enksog approximation, Eq. (E21). The time correlation function 
is (see Eq. (E15)) 

fig 7 ( s ) = V- 1 J dx J c u (x) [1 - V w (x)] e-< AE+ ^) 1<: u (x) (E34) 
with the conjugate flux given by 

1 C v(x) = -(a e + ^)M$(x). (E35) 
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Proceeding as in the case of the shear viscosity, it can be shown that 

e - a (A«+$) 7 ^ {v) = e -s{j+$) 1( . v {v) ) (E36) 



where J is the operator defined in Eq. (89). The time correlation function in Eq. (E34) can 
be expressed as 

fig 7 ( s ) = J dv J c u (v) [1 - V m (x)] e-< J+ ^) 1(U (v) . (E37) 
Therefore, this Euler order transport coefficient is given in the Enskog approximation by 

C U = -\ (1 - a) (p - 1) + J dv J (U (v) [1 - V^\x)] V (U (v) , (E38) 
where V^u (v) is the solution to the integral equation 



(V+f) V c u(v)=[l-VV(x)] 7c, (v). 



(E39) 



d. The bulk viscosity k 

The expression for the bulk viscosity reported in ref. [9] is 

k = lim Vl K H (s) = Q K H (0) + lim f ds'Q K G (s ) , (E40) 

Jo 



where 

O 

and 



h ifi) = ~ J dT tr H (1 - V) e-< c ^)M K (T), (E41) 



= J dTtrH(l-V)e-< c ^)r K (T). (E42) 

The direct flux tr H is the trace of the momentum flux given in Eq. (66), and the adjoint 
density is the same as that for ( u , i.e., 

M K (T) = M c u(F) = -£> • -foPhF)- (E43) 

r=l 

Consequently, the conjugate flux is also the same, T K = T^u, given by (E29). In the Enskog 
approximation, the time correlation function fi^ (s) becomes 

n«j(s) ~ V- 1 J dxJ K (x) [1 -V {1 \x)] e-< AE+ ^M^(x), (E44) 
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with the direct flux given by 

v\ , 1 ding® (a) 



£ + 



d? 2d dnh 

(i + ^i 2) W) 

2d 2 



{p-l) 



J dx 2 5 (q 12 - a) 6 (-q 12 ■ gi 2 ) |<Zi2 '912? (v 2 ) (E45) 



y (re), given in Eq. (E32). The instantaneous contribution fi^-(O) is simply 
related to Q V H (0) in (93) through [9] 



n&(o) = ^pn£(o). 



(E46) 



Furthermore, the correlation function fi£(s) vanishes. This can be seen as follows. The 
conjugate flux and the generator of dynamics can be simplified as in the case of ( u above, 
so that they become independent of the coordinate q. Then, the direct flux above can be 
simplified to give 



Jk («i) 



1 + 



[l + a)a d 7r d / 2 g^\a) 



2T (I) d 



d 2 



dn h 



p — 1 

d ' 



(E47) 



Since J K (v) is orthogonal to the subspace spanned by 1 — V^'(x), ^g( s ) vanishes. Therefore, 
the bulk viscosity in the Enskog approximation is simply 



K 



—j- ^(°) = 2d 2 r( d ±±) — ' 1 ' 2 |vi ~ V2 ^ h ^ fh ^ 

Tr-^r (d/2) a(p - 1) 



r(*±!)d 



dv 2 |«i ~v 2 \fj:p {vi)f y h L> (v 2 ) 



,(i) 



(E48) 



e. T/ie thermal conductivity A 



The thermal conductivity is expressed in ref. [9] as 

A = lim Vl^ (s) = n x H (0) + lim f ds Q X G (s') , 

Jo 



where 



and 



Q x H (s 



-{Vd)~ l J dT S(T) ■ (1 - V) e-< c '^Mx(T) 
( s ) = -{Vd)- 1 J dYS{T) ■ (l-V)e'< Z '^)r x (T). 



(E49) 

(E50) 
(E51) 
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The direct flux S is the heat flux in the microscopic balance equation for the energy density, 

N N N 

S = J2 V r V r + YsYs S2 ^' X ^ ' ( E52 ) 



r=l s^r 



with 



s 2 (x r , x s ) = - + - 5(q rs - a)<d (-q ra ■ g rs ) (q rs • 9rs) 2 (q rs • G rs ) q rs , 



(E53) 



where G rs = (r + v s )/2. The conjugate density in Eq. (E50) is 



1 N d 
M x {Y) = --Y j qr^-[v rPh {Y)] 



r=l 



(E54) 



(E55) 



and the associated flux is 

T A (r) = - 

In the Enskog approximation, the time correlation function determining the Helfand form 
of the thermal conductivity A is 



£-|}M A (r). 



where 



Sl$ I (s)~V- 1 J dxJ x (x)- [l-V w (x)] e-< AE -^)M^\x), 



J\\X\) = — 3 j / dx 2 5{q 12 - a)Q(-q 12 - g 12 ) 



(E56) 



d d 
x(qi2 ■ 0i2) 2 (<Zi2 • G 12 )qi2fh 1] (v 2 ) 



and 



M 1} (*i) = 



<?i <9 



(E57) 



(E58) 



2 9wj 

The function J x defined in Eq. (E57) differs from the generic form given in (E10) by a 
velocity independent term which does not contribute to VLjj. The contribution from the 
initial correlation function to the Green-Kubo form of this transport coefficient is given by 

■d-> 

2vr i/2 r ~p±d)' 



^(0) = V ^L 4-1) / dv, j dv 2 ti\ Vl )tf\v 2 ) 



X 



1 3 

^12^12 + #12(012 • Gi 2 f + -g\ 2 + - (212 • G 12 ) 



while the time- dependent correlation function becomes 

3 (p-l)i r .,2 



^ (a) = 



1 + 



d + 2 



Jdv V -^.[l-V^(v)] e-< J -^ lx (v). 



(E59) 



(E60) 
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Proceeding as in the case of the shear viscosity in Sec. VII, it is found that 



7a 



A 



2 J 2 q d^ ' WW ■ 



Thus the thermal conductivity in the Enskog approximation is given by 

3 (p-l)1 r .,2. 



A = ^(0) 



1 + 



d + 2 



where A is the solution to the integral equation 

Co' 



j-2L)ALv)=[l-VV(v)]>y x (v), 



with J the operator defined in Eq. (89). 



/. The coefficient fi 



(E61) 



/'2 
dv^--[l-V {1) (v)]A(v), (E62) 



(E63) 



This is a new transport mechanism for granular fluids that arises due to the inelasticity 
of collisions. As discussed in detail in ref. [9], this transport coefficient consists of two 
time correlation functions, one of which can be recognized as the time correlation part of 
the thermal conductivity. Therefore, the quantity in terms of which the time correlation 
function takes the simplest form is the linear combination 

dln( h 

a = iimu^ ys) = m'h(U) -I- nm 

'o 



d In rih 



X = lim W H (s) = ft£(0) + lim f ds' fi£ (s') , (E64) 

Jo 



with 



and 



K (s) = -(Vd)- 1 / dT S(T) -(1-V) e'^M^r), 



(«) = -(Yd)- 1 / dT S(T) -(1-V) e" s£ Y M (r). 



(E65) 
(E66) 



The direct flux S above is the heat flux given in Eq. (E52). The conjugate density is 



Ai-p = J drr 



dr r 



5n(r) d Inn 5T (r) 



Spih\F\y] \ 
5n{r) 7 Co 



y=Vh 



(E67) 



y=yh 



where the functional derivative with respect to the density in the last equality is to be taken 
at constant cooling rate, as indicated. In the Enskog approximation, the Helfand form 
becomes 



%(s)^V- 1 J dxJ v {x)- [1-V (1 \x)} e-< AB -%)M$\x). 



(E68) 
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Here, J-p(x) = J\ (x), given in Eq. (E57), and 



d In Co d 



(E69) 



9 In nh dv 

The Green-Kubo form is determined from fi#(0) and fi^ (s). Direct evaluation shows the 
relationship 



!%(o) = -2 



mnC ° n» H (o). 



91nn/j 

The time correlation function (s) can be simplified as in the previous cases yielding 



(E70) 



W G (s) =JdvJ 7T (v) ■ [1 - P« («)] e"^ (») , 



(E71) 



with 



i+ l91n^V) 
2 dlnn^ 



Qi [if («) 



o gig Co 
am n/j 



(E72) 



The Green-Kubo form for the transport coefficient fj, is 
alnC. 



/,-2^A = -2 mnCo 



(9 In n/j 9 In 

where C (v) is a solution to the integral equation 



(0) + J dvJ x (v) ■ [1 - V m (x)} C (v) , (E73) 



JC(v)=[l-VW{x)]J x {v). 



(E74) 



All of the above results agree in detail (for d — 3) with those obtained in ref. [11] by 
means of the Chapman-Enskog solution of the non-linear Enskog kinetic equation. 
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